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Abstract 

Due  to  its  heat  integration  utility,  heat  exchange  reformer  (HER)  is  suitable  to  use  for  PEMFC-based  residential  power  generation  system.  Since 
dynamics  response  of  reformer  affects  overall  dynamics  of  fuel  processing  system,  response  of  HER  with  step  change  in  various  input  parameters 
also  needs  to  be  studied.  In  this  study,  we  present  an  improved  distributed  dynamic  model  for  HER  which  is  capable  of  computing  load  following 
characteristics.  The  dynamic  model  consists  of  20  partial  differential  equations  which  are  discritized  using  spline  collocation  on  finite  elements. 
Resulting  differential  equations  along  with  boundary  conditions  are  solved  by  stiff  solver  which  utilizes  variable  order  method  based  on  numerical 
differentiation  formulae.  Step  variations  in  various  input  parameters  are  considered  and  dynamic  response  of  HER  is  simulated  at  those  variations. 
The  results  provide  valuable  information  about  dynamic  nature  of  HER  which  can  be  used  to  frame  suitable  control  strategy. 

©  2006  Published  by  Elsevier  B.V. 
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1.  Introduction 

In  the  last  decade,  interest  in  cleaner  energy  technologies 
such  as  fuel  cells  has  increased  significantly,  particularly  for 
power  generation  and  automobile  applications.  The  fuel  pro¬ 
cessing  system  for  residential  power  generation  system  based 
on  PEMFC  is  explained  in  Fig.  1.  The  incoming  natural  gas 
is  heated  and  desulphurised  first.  The  steam  is  then  added  and 
methane-steam  mixture  is  sent  to  the  reformer.  The  reformate 
gas  is  cooled  and  passes  through  the  first  shift  converter  where 
majority  of  carbon  monoxide  is  converted  to  carbon  dioxide. 
The  remaining  carbon  monoxide  is  removed  in  second  shift 
converter.  To  obtain  allowable  levels  of  carbon  monoxide  the 
hydrogen  rich  gas  is  sent  to  PROX  reactor  where  it  is  oxidized 
to  carbon  dioxide.  The  exit  gas  is  sent  to  fuel  cells  to  generate 
electricity. 

The  high  endothermicity  of  the  reforming  reaction  and  lim¬ 
ited  wall  heat  transfer  coefficients  necessitate  considerable 
energy  input  requirement  in  reformer  [1].  Hence,  in  such  cir¬ 
cumstances,  the  integration  of  heat  duties  plays  an  important 
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role  to  keep  the  reactor  compact  and,  at  the  same  time,  producing 
the  required  amount  of  hydrogen.  The  methane  heat  exchange 
reformer  (HER)  designed  by  Haldor-Topsoe  utilizes  the  concept 
of  integrating  different  heat  duties  which  results  into  complex 
flow  arrangement  for  such  reactor.  Here,  a  combination  of  co¬ 
current  and  countercurrent  flow  arrangements  between  flue  gas 
and  process  gas  is  devised  to  maximize  the  heat  recovery.  The 
HER  has  been  used  before  for  250  kW  Unocal  MCFC  demon¬ 
stration  project  in  USA  and  also  for  Dutch  50  kW  MCFC  system 
[2].  Due  to  its  heat  integration  utility,  such  reactor  can  also 
be  used  to  PEMFC-based  residential  power  generation  system. 
Another  important  requirement  for  reformer  in  the  fuel  process¬ 
ing  system  for  residential  power  generation  system  is  that  they 
should  operate  in  such  a  manner  that  it  shows  rapid  response  to 
the  load  change  and  also  have  flexible  controllability.  This  moti¬ 
vates  to  design  HER  suitable  for  use  in  PEMFC-based  residential 
power  generation  system  and  analyze  its  dynamic  performance. 

A  previous  simulation  study  on  methane  heat  exchange 
reformer  for  MCFC  power  generation  system  simplified  the 
model  by  considering  each  catalyst  bed  as  a  well-mixed 
reactor  [2].  Since  such  simplifications  do  not  represent  the 
fixed  bed  reactors  with  reasonable  accuracy,  we  present  an 
improved  dynamic  model  without  the  assumption  of  the  well- 
mixed  reactor  suitable  for  use  for  PEMFC-based  residential 
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Nomenclature 

cp 

heat  capacity  (Jkg-1  K-1) 

cPi 

heat  capacity  of  component  i  (Jkg-1  K-1) 

C 

concentration  (molm-3) 

d 

hydraulic  diameter  (m) 

EP 

activation  energy  (J  mol-1) 

h 

wall  heat  transfer  coefficient  (W  m-2  K-1) 

Hi 

enthalpy  of  ith  species  (J  mol-1) 

AH, 

heat  of  reaction  of  yth  reaction  (J  mol-1) 

A//298 

heat  of  reaction  at  298  K  (J  mol-1) 

h 

thermal  conductivity  of  gas 

kw 

thermal  conductivity  of  wall 

k 

rate  constant  of  reaction  i 

Kej 

equilibrium  constant  of  yth  reaction 

L 

length  of  each  catalyst  bed  (m) 

Nu 

Nusselt  number  (=dh/kg) 

P 

partial  pressure  (Pa) 

R 

universal  gas  constant  (J  mol-1  K-1) 

Re 

Reynolds  number  ( -dupjkg ) 

P\.i 

rate  of  reaction  i  (mol  s-1  kg-1  catalyst) 

rr,i 

rate  of  generation  of  component  i  in  reformer 
(mols-1  kg-1  catalyst) 

Sr 

heat  transfer  area  per  unit  volume  (m2/m3) 

t 

time  (s) 

thickness  of  wall  (m) 

T 

temperature  (K) 

Pref 

reference  temperature  (298  K) 

U 

superficial  velocity  (ms-1) 

z 

length  of  reactor  (m) 

Greek  symbol 

P 

density  (kg  m  3) 

Subscripts 

g 

gas 

i 

ith  reactant 

j 

y'th  reaction 

f 

flue  gas 

fl 

first  flue  gas  channel 

f2 

second  flue  gas  channel 

pi 

process  gas  channel 

p2 

product  gas  channel 

r 

catalyst  bed 

rl 

first  catalyst  bed 

r2 

second  catalyst  bed 

w 

wall 

wl 

first  wall 

w2 

second  wall 

w3 

third  wall 

w4 

fourth  wall 

solved  numerically  and  dynamic  performance  of  HER  for  step 
variations  in  various  input  parameters  is  simulated. 

2.  Reactor  geometry 

As  shown  in  Fig.  2  [3],  the  Haldor-Topsoe  HER  comprises 
of  pressure  vessel,  two  catalyst  bed  and  burner.  The  reformer 
fuel  is  combusted  in  the  burner  located  at  the  bottom  of  the 
pressure  vessel.  Feed  gas  flows  downwards  through  first  catalyst 
bed  where  it  is  heated  by  convection  from  both  the  flue  gas  and 
the  reformed  product  gas,  both  flowing  countercurrent  to  the 
feed.  After  leaving  the  first  catalyst  bed,  the  partially  reformed 
gas  is  fed  to  the  top  of  second  catalyst  bed  where  it  is  heated 
by  flue  gas  and  the  partially  reformed  process  gas  as  shown  in 
Fig.  2.  Both  catalyst  beds  are  filled  with  Ni/Mg0-Al203  catalyst 
where  steam  reforming  of  methane  takes  place.  It  is  also  clear 
from  Fig.  1  that  the  section  of  the  reactor  to  the  left  of  the  burner 
channel  consists  of  six  channels  including  two  channels  for  fixed 
bed  reactor,  two  channels  for  the  process  gas  and  the  outer  two 
channels  for  the  flue  gas  flow.  Each  channel  of  the  reactor  is 
considered  to  have  rectangular  cross  section.  Due  to  symmetry 
of  the  reactor,  only  one  section  of  the  reactor  is  modeled. 

3.  Reaction-scheme  and  kinetics 

3.1.  Reformer 

Steam  reforming  of  methane  takes  place  in  presence  of  the 
Ni/Mg0-Al203  catalyst  which  is  described  by  following  reac¬ 
tion. 

CH4  +  H2O  **  CO  +  3H2,  A#298  =  206 kJ  mol-1  (1) 

CO  +  H20  CO2  +  H2,  AH 2 98=  —41  kJ  mol-1  (2) 

CH4  +  2H20  **  C02  +  4H2,  A// 298  =  165kJ mol-1  (3) 

Kinetic  information  for  these  reactions  is  taken  from  Xu  and 
Froment  [4]. 


^r.  l  = 


k2 


Pcn4  Pr2o  ~ 


ph2  pco  \  1 


K, 


ft r.2  =  -  -  (  ftcoftH2Q  - 


ftr,3  = 


ei 

ftlh  pco 


Ql 


(4) 


1 


K, 


?2 


Ql 


(5) 


p  p2  pn2  pco2 
Pcu4  "h20  “  " 


1 


KQ3  /  £2r 


(6) 


where 

Qr  =  1  +  KqoPco  +  ^h2Ph2  +  Kch4Pch4  + 
3.2.  Burner 


^h2o^h2o 

P\h 


power  generation  system.  The  model  consists  of  12  partial 
differential  equations  which  are  descritized  using  orthogonal 
spline  collocation  on  finite  elements.  The  resulting  ordinary 
differential  equations  along  with  boundary  conditions  are 


Reaction  scheme  for  the  burner  is  as  below 
CH4  +  2O2**  C02  +2H20,  AH 29s  =  -802.3  kJmor1 


(7) 
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Fig.  1.  Fuel  processing  system. 


0.5C>2  +  H2  H2O,  AH298  =  — 241.8  kJ mol  *  (8) 

CO  +  0.502  **  CO2,  AH29&  =  -283  kJ  mol-1  (9) 

The  rate  equations  for  catalytic  combustion  are  not  included 
in  the  model  equations.  Assuming  that  combustion  reactions 
are  comparatively  faster  than  stem  reforming  reaction,  flue  gas 
temperature  and  composition  at  burner  outlet  can  be  obtained 
from  equilibrium  calculations. 


First  catalyst  bed: 


3Cri  J  3Cr]  i 

Wf  —  ~Uv  o  ^  Prrrl,h 
ot  OZ 

Second  catalyst  bed: 

dCr2,i  dCv2J  . 

—  ~aY  —  b  Pvrv2,ii 


=  1-5 


dt 


3  z 


=  1-5 


4.2.  Energy  balance  equations 


(10) 


(11) 


4.  Simulation  model 

Both  fixed  bed  catalyst  beds  are  modeled  as  one-dimensional 
pseudo-homogenous  reactor  of  the  equal  length  and  of  equal 
catalyst  bed  density.  In  order  to  simplify  the  model,  following 
assumptions  are  made: 

•  Pressure  inside  both  the  flue  gas  and  the  reformer  channels  is 
assumed  to  be  constant. 

•  Conduction  and  convection  are  assumed  to  be  predominant 
heat  transfer  mechanism  while  radiation  effect  is  neglected. 

•  Heat  loss  to  surroundings  is  neglected. 

•  Both  the  process  gas  and  the  fuel  gas  are  treated  as  ideal  gas. 

4.1.  Mass  balance  equations 


Similarly,  energy  balance  equations  are  given  by  partial  dif¬ 
ferential  equations.  Here,  the  first  and  the  second  wall  refer  to 
the  two  opposite  sides  of  the  second  catalyst  bed  along  which 
the  gas  flow  occurs.  Similarly,  the  third  and  the  fourth  wall  form 
the  opposite  sides  of  the  second  catalyst  bed.  The  first  flue  gas 
channel  and  the  first  process  gas  channel  are  on  the  opposite 
sides  of  the  second  catalyst  bed  while  the  second  flue  gas  chan¬ 
nel  and  the  second  preheat  channel  are  on  the  opposite  sides  of 
the  first  catalyst  bed  as  shown  in  Fig.  3.  With  this  introduction  of 
the  nomenclature,  the  energy  balance  equations  can  be  written 
as  follows: 

•  First  catalyst  bed: 

3Trj  3Tri 

^vCprl~di~  ~  ~u^>vCP^~^z  Ay  v^rl,j(— A/frij) 

7=1 

+  hY\SY\(Tw4  —  TY  i)  +  ^rl5rl(Tw3  —  TY\) 


Based  on  these  assumptions  stated  before,  mass  balance  equa¬ 
tions  can  be  expressed  as  partial  differential  equations  as  follows: 


(12) 


Second  catalyst  bed: 


371-2  371-2  v-^ 

Pvc pr2  —  —  upYCpr2— —  b  Pr/  jRy2, j(~ A77r2, j) 

7=1 


3 1  3  z 

+  hY2SY2(Tw\  -  TY2)  +  hY2 5r2(7w2  -  71-2) 


(13) 


Flue  gas  channel  1 : 

37fi  37fi 

PfcPf~Z  =  UfPfCp^— —  b  AfiSfi(rwi  Tfi ) 

ot  oz 

Flue  gas  channel  2: 

37f2  37f2 

PfOpf~Z  =  UfPfCp^— -  b  hf2 Sf2(T'w4-  Tfe) 

ot  oz 

Process  gas  channel: 

3Tpi  3Tpi 

Pp\cpvi  ~  ~ u*PpicPpi  ~ 1“  ^pi*Spi(7w2  —  7"pi) 


(14) 


(15) 


Fig.  2.  Schematic  diagram  of  the  reactor  [3]. 


(16) 
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Fig.  3.  Dynamic  behavior  of  HER  for  15%  step  increase  in  inlet  process  gas  flow:  (a)  product  gas  composition,  (b)  product  gas  flow  rate  and  (c)  exit  gas  temperatures. 


•  Product  gas  channel: 


•  Fourth  wall: 


dTp2 

Pplcpp2  dt 


3rp2 

MrPp2Cpp2  ^ 


T  //p2.Sp2^w3 


97w4 

Pwcpv •  -  — k 


dt 


w 


d2Tw4  hf2  /  .  hv\ 

w4  +  — (rf2-rw4)  +  — (rrl 


dz2 


t 


w 


t 


w 


Tw  4) 

(21) 


First  wall: 

3Fwi  d2Tvi\  hfi  hr  1 

PwCpw  f—=kv/  +  —  (7fl  -  rwi)  +  —  (Tri  -  rwl) 


dt  W  dz 2  f 


w 


t 


w 


(18) 


Second  wall: 

97*2  9“  7*2  /ipl  /?rl 

PwCpw— ?=*w^-F+— (rpi  -  7*2)  +  —  (Til  -  rw2) 
of  dz  f\v 

(19) 

Third  wall: 


97’w3  92Tw3  /zP2  /tr2 

PwCpw— ^=^w^r^  +  -^(7p2  -  7^3)  +  -f{Tt2  -  rw3) 

07  oz  r\v  r\v 

(20) 


The  following  boundary  conditions  are  added  to  the  modeling 
equations. 

Here,  the  top  of  each  channel  and  each  catalyst  bed  is  con¬ 
sidered  as  z  =  0  and  the  bottom  as  z  =  L. 

At  z  =  0 

dTw\(t)  _  dTw2(t)  _  dTwi(t)  _  dTw4(t) 
dz  dz  dz  dz 

Quit)  =  CrU(  0),  i  =  1-5 

Cr2,z(Olz=0  =  Cr\j(t)\z=L,  i  =  1-5 

TMt)  =  Tr\  (0) 


(22) 

(23) 

(24) 

(25) 


7r2(Olz=0  =  Tp\(t)\z=L 


(26) 
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Boundary  conditions  at  z  =  L 

3Twi(Q  _  d?w2 (0  _  dTw3(t)  _  dTw4(t) 
3  z  dz  3  z  3  z 


integrator  that  utilizes  variable  order  solver  based  on  numeri¬ 
cal  differentiation  formulae.  The  computations  are  performed  in 
(27)  MATLAB  programming  environment. 


Tr2(t)\z=L  =  Tp2(t)\z=L 


(28) 


5.  Results  and  discussion 


4.3.  Heat  transfer  coefficients 


The  reformer  wall  heat  transfer  coefficient  is  estimated  from 
correlation  (29)  [5]. 

Nu  =  0.027 Cc(Re)0S(Prf33  (29) 


where  Cc  is  constant  which  depends  on  size  of  the  catalyst  par¬ 
ticles. 

Depending  on  the  type  of  the  flue  gas  flow,  the  heat  transfer 
coefficient  can  be  estimated  from  following  correlations  [6]. 
For  laminar  flow  of  the  flue  gas: 


f  df'3 

Nu  =  1.86  [RePr-  ) 


(30) 


For  turbulent  flow  of  the  flue  gas: 

Nu  =  0.027 (Re)0S{Prf33  (31) 


The  base  case  reformer  configuration  is  shown  in  Table  1. 
Simulation  results  for  this  configuration  shows  that  HER  gen¬ 
erates  4.32  mol  min-1  of  hydrogen  flow  which  corresponds  to 
6.96  kW  power  assuming  40%  efficiency  for  fuel  cell  system. 
This  is  taken  as  reference  state  and  following  step  variations  are 
considered  to  study  dynamic  performance  of  HER: 

•  ±15%  step  changes  in  inlet  process  gas  flow  rate; 

•  ±15%  step  changes  in  process  gas  pressure; 

•  ±15%  step  changes  in  inlet  burner  flue  gas  pressure  and  flue 
gas  flow  rate; 

•  ±15%  step  changes  in  steam  to  methane  ratio; 

•  step  change  in  process  gas  inlet  temperature  from  400  to 
460  K. 

5.1.  Dynamics  of  process  gas  flow  rate 


4.4.  Fluid,  catalyst  and  wall  properties 

Thermal  conductivity  of  fuel  gas  can  be  expressed  as  a  func¬ 
tion  of  temperature  from  (32)  [7]. 

kg  =  1.679  x  10"2  +  5.073  x  10_57’  (32) 

Enthalpy  of  component  /in  reformer  gas  is  given  by  (33): 

Hi  =  f  cPi  d T  (33) 

2  Tef 

Heat  capacity  of  gas  species  is  evaluated  as  a  function  of  tem¬ 
perature  from  (34): 

c pi  =  a  ±  bT  ±  cT2  (34) 

Viscosity  of  reformer  gas  mixture  is  evaluated  as  function  of  con¬ 
centration  and  temperature  using  corresponding  states  methods 
[8]. 

Thermal  conductivities  and  heat  capacities  of  reformer  cat¬ 
alyst  is  taken  as  that  of  alumina  [6].  Thermal  conductivity  and 
heat  capacity  of  wall  is  taken  as  that  of  iron  [6].  Thermal  con¬ 
ductivities  of  individual  gas  species  is  taken  from  [9].  From 
component  thermal  conductivity,  mixture  thermal  conductivity 
can  be  estimated  using  method  of  Linsay  and  Bromley  [10]. 

4.5.  Numerical  technique 

The  set  of  20  partial  differential  equations  (( 10)— (21))  are  dis- 
critized  using  orthogonal  spline  collocation  on  finite  elements 
[11].  The  reactor  length  is  divided  into  10  intervals  with  two 
collocation  points  in  each  interval.  Discritized  partial  differen¬ 
tial  equations  along  with  boundary  conditions  form  system  of 
stiff  differential  algebraic  equations  which  are  solved  using  stiff 


The  residential  fuel  cell  power  system  is  often  required  to 
operate  under  partial  load.  To  enable  to  get  different  power  out¬ 
put,  manipulation  of  inlet  process  gas  flow  of  the  reformer  can 
be  one  of  the  possible  variations.  Hence,  dynamic  performance 
of  the  HER  under  ±15%  step  variations  in  inlet  process  gas  flow 
is  analyzed  here. 

Simulated  composition  profiles  of  the  product  gas  for  step 
increase  in  process  gas  flow  rate  are  shown  in  Fig.  3(a).  Results 
indicate  that  methane  and  steam  compositions  increase  while 
hydrogen  composition  decreases  with  time.  This  trend  is  similar 
to  the  experimentally  observed  dynamic  concentration  profile 
[12].  The  hydrogen  composition  decreases  from  49.72%  to 
45.35%  with  net  gain  of  —4.37%.  Similarly,  methane  concentra¬ 
tion  increases  from  3.85%  to  6.1 1%  with  net  gain  of  2.26%.  For 
step  decrease  in  process  gas  flow  rate,  it  is  found  that  methane 
composition  decreases  from  3.85%  to  1.74%  with  net  gain  of 
—2.1 1%  while  hydrogen  composition  increases  from  49.72%  to 
53.3%  with  net  gain  of  3.58%.  Results  show  that  all  responses 

Table  1 


Design  specifications  and  operating  conditions  for  the  base  case 


Design  specification/operating  condition 

Value 

Length  of  each  catalyst  bed 

1  m 

Reformer  channel  thickness 

0.01  m 

Thickness  of  each  flue  gas  channel  thickness 

0.01  m 

Thickness  of  each  preheat  channel  thickness 

0.01  m 

Reformer  catalyst  bed  density 

100  kg  m~3 

Wall  thickness 

1  x  10-3  m 

Reformer  pressure 

2  bar 

Burner  pressure 

8  bar 

Steam/methane  ratio 

3 

Reformer  inlet  gas  temperature 

400  K 

Reformer  inlet  gas  flow 

0.1083  mols-1 

Flue  gas  flow 

0.6873  mol  s_1 

508 


K.S.  Patel,  A.K.  Sunol  /  Journal  of  Power  Sources  161  (2006)  503-512 


are  stable  without  any  overshoot  and  reaches  stead  state  within 
500  s. 

For  step  increase  in  process  gas  flow  rate,  the  total  molar 
flow  rate  increases  from  4.34  to  4.79  mol  min-1  with  net  gain 
of  10.36%  and  hydrogen  flow  rate  increases  from  2.16  to 
2.174  mol  min-1  with  net  gain  of  0.648%  as  shown  in  Fig.  3(b). 
Similarly,  the  step  decrease  in  flow  rate  decreases  in  total 
molar  flow  rate  to  3.838  mol  min-1  and  hydrogen  flow  rate 
to  2.046  mol  min-1  with  net  gains  of  —11.65%  and  —5.28%, 
respectively.  Also,  its  should  be  noted  that  for  both  cases,  flow 
rates  increases  rapidly,  reaches  maximum  and  then  decreases 
with  slight  overshoot.  This  dynamic  behavior  of  HER  matches 
in  quality  with  the  previous  simulation  studies  [2] .  Higher  gain 
in  total  molar  flow  rate  suggests  that  the  step  change  in  inlet 
flow  rate  is  suitable  strategy  for  manipulating  total  flow  rate  of 
product  gas. 

Fig.  3(c)  shows  thermal  behavior  of  HER.  For  step  increase 
in  flow  rate,  product  gas  temperature  increases  initially  from  826 
to  836.6  K  and  then  decreases  to  steady  state  value  of  806.9  K 
with  net  gain  of  —2.31%.  Flue  gas  temperature  decreases  from 
1066  to  1038  K  with  net  gain  of  —2.63%.  Similarly,  on  step 
decrease  in  process  gas  flow,  the  process  gas  temperature  first 


decreases  to  8 13.6  K  and  then  increases  to  845.7  K  with  net  gain 
of  2.39%.  Flue  gas  temperature  increases  to  1098  K  with  net  gain 
of  3%.  Although,  all  the  temperature  responses  reach  steady  state 
within  500s,  overshoots  in  thermal  response  can  affect  overall 
performance  of  fuel  processing  system. 

5.2.  Dynamics  of  process  gas  pressure 

Pressure  change  is  one  of  the  important  parameters  that  can 
affect  desirable  operation  of  HER.  Therefore,  dynamic  perfor¬ 
mance  of  HER  with  ±15%  step  increase/decrease  of  pressure  is 
analyzed  in  this  section. 

As  shown  in  Fig.  4(a),  15%  step  decrease  in  process  gas  pres¬ 
sure  decreases  methane  molar  composition  from  3.85%  to  1 .49% 
with  net  gain  of  —2.36%  while  hydrogen  composition  increases 
from  49.72%  to  53.8%  with  net  gain  of  4.08%.  Similarly,  for 
step  increase  in  process  gas  pressure,  methane  concentration  is 
found  to  increase  from  3.85%  to  6.33%  with  net  gain  of  2.48% 
and  hydrogen  composition  decreases  from  49.72%  to  44.85% 
with  net  gain  of  —4.87%.  Results  show  that  all  responses  are 
stable  without  any  overshoot  and  reaches  stead  state  within 
500  s. 
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Fig.  4.  Dynamic  behavior  of  HER  for  15%  step  decrease  in  process  gas  pressure:  (a)  product  gas  composition,  (b)  product  gas  flow  rate  and  (c)  exit  gas  temperatures. 
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For  step  decrease  in  process  gas  pressure,  the  total  molar 
flow  rate  decreases  from  4.34  to  3.86  mol  min-1  with  net  gain 
of  —11.06%  and  hydrogen  flow  rate  decreases  from  2.16  to 
2.076  mol  min-1  with  net  gain  of —3.89%  as  shown  in  Fig.  4(b). 
Similarly,  the  step  increase  in  process  gas  pressure  causes 
increase  in  total  molar  flow  rate  to  4.78  mol  min-1  with  net  gain 
of  10. 14%.  The  hydrogen  flow  rate  decreases  to  2. 14  mol  min-1 
with  net  gain  of  —0.93%.  Also,  is  should  be  noted  that 
for  both  cases,  flow  rates  changes  rapidly,  reaches  extreme 
value  and  then  approaches  the  steady  state  value  with  slight 
overshoot. 

Fig.  4(c)  shows  thermal  behavior  of  HER.  For  step  decrease  in 
process  gas  pressure,  product  gas  temperature  decreases  initially 
from  826  to  813.6  K  and  then  increases  to  steady  state  value  of 
845.5  K  with  net  gain  of  2.39%.  Flue  gas  temperature  increases 
from  1066  to  1098  K  with  net  gain  of  3%.  Similarly,  on  step 
increase  in  process  gas  pressure,  the  process  gas  temperature 
first  increases  to  836.4  K  and  then  decreases  to  810.6  K  with  net 
gain  of  — 1.86%.  Flue  gas  temperature  decreases  to  1038  K  with 
net  gain  of  —2.63%.  It  is  found  that  all  the  temperature  responses 
reach  steady  state  within  500s. 


5.3.  Dynamics  of  inlet  process  gas  temperature 

Since  inlet  process  gas  is  heated  and  partially  reformed  in  first 
bed  in  HER,  inlet  temperature  plays  important  role  in  overall 
performance  of  HER.  Hence,  dynamics  of  step  increase  of  inlet 
process  gas  temperature  from  400  to  460  K  is  studies  in  this 
section. 

As  shown  in  Fig.  5(a),  product  gas  composition  dynamics 
with  step  increase  in  process  gas  temperature  is  faster  that  earlier 
discussed  variations.  Hydrogen  concentration  increases  from 
49.72%  to  55.54%  and  then  approaches  to  steady  state  value 
of  51.07%  with  net  gain  of  1.35%.  Similarly,  methane  concen¬ 
tration  decreases  from  3.85%  to  0.269%  and  approaches  steady 
state  value  of  3. 1 1  %  with  net  gain  of  —0.74%.  Results  also  reveal 
that  the  composition  response  is  faster  and  reaches  steady  state 
in  about  10  s. 

Fig.  5(b)  explains  the  flow  rate  variations.  The  total  molar 
flow  rate  decreases  from  4.34  to  3.78  mol  min-1,  increases  to 
4.365  mol  min-1  and  then  decreases  to  3.82  mol  min-1  with 
net  gain  of  —11.98%.  Hydrogen  molar  flow  rate  shows  similar 
behavior.  It  decreases  from  2.16  to  1.88  mol  min-1,  increases 
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Fig.  5.  Dynamic  behavior  of  HER  for  step  increase  in  process  gas  temperature  from  400  to  460  K:  (a)  product  gas  composition,  (b)  product  gas  flow  rate  and  (c)  exit 
gas  temperatures. 
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to  2.4  mol  min-1  and  then  decreases  to  1 .95  mol  min-1  with  net 
gain  of  —9.72%.  Also,  is  should  be  noted  that  for  both  cases, 
flow  rates  changes  rapidly,  with  the  steady  state  value  reaches 
only  in  about  10  s. 

As  shown  in  Fig.  5(c),  product  gas  temperature  increases 
initially  from  826  to  950  K  and  then  decreases  to  steady  state 
value  of  839  K  with  net  gain  of  1.57%.  Flue  gas  temperature 
decreases  from  1066  to  1083  K  with  net  gain  of  1.59%.  It  is 
found  that  thermal  response  of  HER  is  slower  than  flow  response 
and  reaches  steady  state  values  in  about  500  s. 

5.4.  Dynamics  of  flue  gas  pressure  and  flue  gas  flow  rate 

Since  steam  reforming  is  highly  endothermic  reaction,  flue 
gas  pressure  plays  important  role  to  maintain  high  temperature 
and  thus  creates  favorable  condition  for  steam  reforming.  Here, 
dynamic  performance  of  HER  for  step  variation  of  flue  gas  pres¬ 
sure  is  analyzed. 

As  shown  in  Fig.  6(a),  step  increase  in  flue  gas  pressure 
decreases  methane  composition  from  3.85%  to  2.32%  with 
net  gain  of  —1.53%  and  increases  hydrogen  composition  from 


49.72%  to  52.3%  with  net  gain  of  2.58%.  Similarly  for  step 
decrease  in  burner  pressure,  net  gains  in  methane  and  hydrogen 
compositions  are  found  to  be  —4.01%  and  2. 13%.  Results  reveal 
that  all  responses  are  stable  without  any  overshoot  and  reaches 
stead  state  within  500s. 

For  step  increase,  the  total  molar  flow  rate  increases  from 
4.34to  4.46  mol  min-1  with  net  gain  of  2.77%  and  hydrogen 
flow  rate  increases  from  2.16  to  2.33  mol  min-1  with  net  gain 
of  7.87%  as  shown  in  Fig.  6(b).  Same  way,  the  step  decrease 
in  burner  pressure  causes  net  gain  of  —3.45%  and  —11.3% 
in  total  and  hydrogen  molar  flow  rate,  respectively.  Responses 
indicated  that  steady  state  is  reached  in  about  500s.  The  gain 
values  also  suggest  that  fluctuation  in  burner  inlet  pressure 
effects  hydrogen  molar  flow  rate  more  than  the  total  molar  flow 
rate. 

Fig.  6(c)  shows  thermal  behavior  of  system.  For  step  increase 
in  flue  gas  pressure,  product  gas  temperature  increases  from 
826  to  852  K  with  net  gain  of  3.15%.  Flue  gas  temperature 
increases  from  1066  to  1098  K  with  net  gain  of  3%.  Similarly, 
step  decrease  in  burner  pressure  causes  decrease  in  product  gas 
temperature  and  flue  gas  temperature  with  net  gains  of  —3.63% 
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Fig.  6.  Dynamic  behavior  of  HER  for  15%  step  increase  in  flue  gas  pressure:  (a)  product  gas  composition,  (b)  product  gas  flow  rate  and  (c)  exit  gas  temperatures. 
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Fig.  7.  Dynamic  behavior  of  HER  for  15%  step  increase  in  steam  to  methane  ratio:  (a)  product  gas  composition,  (b)  product  gas  flow  rate  and  (c)  exit  gas  temperatures. 


and  —3.66%,  respectively.  All  the  temperature  responses  reach 
steady  state  within  500  s  without  any  overshoot.  Simulated 
results  for  step  variation  in  flue  gas  flow  rate  show  similar  gains 
for  compositions,  flow  rates  and  temperatures. 

5.5.  Dynamics  of  steam  to  methane  ratio 

Dynamic  performance  of  HER  with  step  increase/decrease 
of  steam  to  methane  ratio  is  analyzed  in  this  section. 

For  15%  step  increase  in  steam  to  methane  ratio,  overshoots 
are  observed  in  compositions  of  all  components  as  shown  in 
Fig.  7(a).  Hydrogen  composition  decreases  from  49.72%  to 
48.51%  with  net  gain  of  —1.21%.  Methane  composition  also 
decreases  from  3.85%  to  2.54%  with  net  gain  of  —1.31%.  It  is 
observed  that  the  response  is  fast  and  reaches  steady  state  value 
in  just  5  s.  Similarly,  for  15%  step  decrease  in  steam  to  methane 
ratio,  hydrogen  and  methane  compositions  are  found  to  increase 
to  50.64%  and  5.65%,  respectively. 

As  with  compositions,  overshoots  are  also  observed  in  total 
molar  flow  rate  and  hydrogen  molar  flow  rate  as  shown  in 
Fig.  7(b).  Total  flow  rate  decreases  from  4.34  to  4.303  mol  min-1 
with  net  gain  of  —0.85%  and  hydrogen  molar  flow  rate  decreases 


from  2.16  to  2.088  mol  min-1  with  net  gain  of  —3.33%.  It 
is  found  that  both  the  responses  are  fast  with  major  response 
achieved  only  in  5  s.  For  step  decrease  in  stem  to  methane  ratio, 
total  flow  rate  increases  from  4.34  to  4.38  mol  min-1  with  net 
gain  of  0.69%  and  hydrogen  flow  rate  increases  from  2.16  to 
2.22  mol  min-1  with  net  increase  of  2.6%. 

As  shown  in  Fig.  7(c),  step  increase  in  steam  to  methane  ratio 
causes  increase  in  product  from  826  to  834  K  with  net  gain  of 
0.97%  and  increase  in  flue  gas  temperature  from  1066  to  1070  K 
with  net  gain  of  0.38%.  For  step  decrease  in  steam  to  methane 
ratio,  product  gas  temperature  is  found  to  decrease  to  818.3  K 
with  gain  of  —0.93%  and  flue  gas  temperature  also  decreases  to 
1061  K  with  gain  of  —0.47%. 

6.  Conclusion 

In  this  work,  we  numerically  investigated  dynamic  perfor¬ 
mance  of  methane  heat  exchange  reformer  that  is  suitable  to 
be  used  for  PEMFC-based  residential  power  generation  sys¬ 
tem.  Improved  distributed  model  for  HER  is  presented  and 
dynamic  effects  of  step  variations  in  various  input  variables 
are  simulated.  Model  equations  form  coupled  partial  differen- 
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tial  equations  which  are  descritized  using  spline  collocation  on 
finite  elements  and  solved  numerically.  The  simulation  results 
indicate  higher  gains  in  total  molar  flow  rate  by  step  changes 
in  process  gas  flow  rate  and  process  gas  pressure  with  slight 
overshoot.  Product  gas  temperature  shows  overshoot  for  step 
changes  in  flow  rate  or  pressure  of  process  gas.  Higher  gains 
in  hydrogen  flow  rate  and  comparatively  lower  gains  in  total 
molar  flow  rate  are  obtained  with  variations  in  flue  gas  velocity 
or  flue  gas  pressure.  Furthermore,  product  gas  temperature  does 
not  show  any  overshoot  for  step  changes  in  flue  gas  pressure 
or  flue  gas  velocity.  Response  time  for  steady  states  is  found 
to  be  500  s  for  all  these  variations.  Step  changes  in  inlet  pro¬ 
cess  gas  temperature  show  rapid  (steady  state  value  reaches  in 
10  s)  response  and  higher  gains  in  total  and  hydrogen  flow  rates 
with  overshoots.  The  product  gas  temperature  for  this  varia¬ 
tion  shows  overshoot  and  reaches  steady  state  in  500  s.  Flow 
rates  show  overshoots  for  step  variation  in  steam  to  methane 
ratio.  Also,  product  gas  temperature  and  total  flow  rate  are 
found  to  be  less  sensitive  to  variations  in  steam  to  methane 
ratio. 
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